Point-defect haloing in curved nematic films 
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We investigate the correlation between the point disclination energies and the surface curvature 
modulation of nematic liquid crystal membranes with a Gaussian bump geometry. Due to the 
correlation, disclinations feel an attractive force that confines them to an annulus region, resulting 
in a halo distribution around the top of the bump. The halo formation is a direct consequence of 
the nonzero Gaussian curvature of the bump that affects preferable configurations of liquid crystal 
molecules around the disclination core. 
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I. INTRODUCTION 

It is generally assumed that topological defects in liq- 
uid crystal membranes are excluded in the ground state 
because they increase the elastic free energy of the system 
to more than that in the defect-free equilibrium state. 
The situation is quite different when the liquid crystal 
phases are embedded in a curved surface. In a curved 
membrane with nonzero Gaussian curvature, the appear- 
ance of topological defects can be favorable or even in- 
evitable because of geometrical and/or topological con- 
straints. A typical example is a disclination-containing 
equilibrium state in a nematic liquid crystal that is con- 
strained in a spherical surface [1]. The defect structure 
in such a nematic spherical shell has recently been recon- 
sidered both theoretically [2-5] and experimentally [6, 7], 
and these investigations revealed richer physical proper- 
ties than those that have initially been suggested. In ad- 
dition to the nematic cases, the smectic phases of liquid 
crystal films [8, 9] and columnar phases of a block copoly- 
mer assembly [10] have also exhibited an interplay be- 
tween disclination structures and the geometry /topology 
of the constraint surface. In particular, smectic ordering 
on a sinusoidal substrate has indicated a direct coupling 
between the substrate curvature and the preferential lo- 
cations of disclinations in equilibrium states [9] . The cou- 
pling is purely a geometric curvature effect in the sense 
that the system involves no topological requirement for 
the existence of disclination. The geometric curvature 
also strongly influences the orientational order of the spin 
lattice models denned on curved surfaces [11-18]. 

Many intriguing results have been obtained with re- 
gard to the global lowest-energy state on a curved surface 
with a minimum number of disclinations, which should 
be relevant to liquid crystal phases at very low tempera- 
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tures. In actual systems, however, a slight physical (and 
chemical) disturbance causes excess disclinations [19-22]. 
These disclinations can travel over a surface via thermal 
excitation, and they flow almost freely if the disclination 
density is sufficiently low to ignore their interaction. The 
only possible force that restricts the free motion of the 
excess disclinations arises from a coupling between the 
surface curvature and the disclination energy. The sub- 
strate curvature breaks the translational symmetry of the 
molecular configuration in the vicinity of the disclination 
core, which may cause a shift in the disclination energy. 
As a result, the energetically preferential positions of ex- 
cess disclinations are believed to be tunable by imposing 
an appropriate curvature modulation. 

In this article, we study the energies of disclinations 
in a nematic thin film assigned on a Gaussian bump. It 
is found that disclinations with a +1 charge prefer posi- 
tioning at inflection points of the bump, implying the for- 
mation of a "halo," i.e., an annulus distribution of excess 
disclinations around the top of the bump. In addition, 
disclinations with a — 1 charge are found to exhibit double 
concentric halo structures, provided the bump height is 
relatively large. These halo formations are consequences 
of the non-trivial coupling between the disclination ener- 
gies and the surface curvature modulation. 



II. MODEL 

We consider nematic liquid crystals confined in a 
thin curved membrane with thickness h, in which the 
molecules align parallel to the tangential directions inde- 
pendently of the depth in the membrane. In line with 
the continuum theory [19, 23], the elastic free energy of 
the molecular ensemble within the volume Ah is given by 

F^hJjA . n ) 2 + ^(Dxn) 2 }, (1) 
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where the unit vector n (called the director field) de- 
scribes the local molecular alignment. D is a vector op- 
erator whose components are the covariant derivative 
[24] on the curved surface. K\ and if 3 are the elastic con- 
stants associated with splay and bend distortions of the 
molecular configuration, respectively [25]. 

Suppose that the director lies within a thin curved shell 
having a Gaussian bump geometry. Points on the bump 
are represented by a three-dimensional vector R(r, (j)) 
given by 

f r 2 \ 

R(r, (j)) = r cos <j)i + r sin <j)j + v exp I — \ fc, (2) 

where i : j, and k are basis vectors of the Cartesian co- 
ordinate system, and r and <\> are the polar coordinates 
assigned on the underlying flat plane (see Fig.l). The ge- 
ometric deviation of a bump from a flat plane is charac- 
terized by a dimensionless parameter a = v/ro that rep- 
resents the ratio of the height to the width of the bump. 
Note that the constant ro in the last term in Eq.(2) de- 
termines the points of inflection of the Gaussian bump, 
i.e., across which the sign of (d 2 R/dr 2 ) • k changes. We 
demonstrate below that ro plays a prominent role in de- 
scribing the correlation between the disclination energy 
and the substrate curvature. 

Let us locate a disclination core on a point R = 
R{r Cl C ), where <j) c is an arbitrary constant without loss 
of generality. The disclination size is sufficiently smaller 
than the bump width, and the perimeter of the discli- 
nation is defined by the locus of points whose distances 
from the core are equal in the sense of geodesic distance; 
therefore, the perimeter slightly deviates from a circular 




FIG. 1: Schematic of a Gaussian bump and mathematical no- 
tations used in the paper. The thick arc depicted on the bump 
is a projection of the thick line segment on the underlying flat 
plane. 



shape when it is placed on a curved surface, whereas it 
is an exact circle when placed on a flat plane. 

For later use, we introduce a new polar coordinate sys- 
tem (/, 0) on the flat plane such that the coordinate origin 
(i.e., I = 0) locates at the point of interaction of the flat 
plane and the vertical line penetrating the disclination 
core (see Fig.l). Then, the basis vectors a± = dR/dl 
and a 2 = dR/dO span the curved surface and allow us 
to write n = n l ai (i = 1,2), where n l and ai depend 
on I and (See Appendix B). The director orientation 
around the core is given by the angle ip between n and 
a\. It is natural to assume that 

fl> = (s-l)0 + fh, (3) 

where s denotes the disclination charge that quantifies 
how many times n rotates about the core. When s = 
+1, for instance, all n contained within the disclination 
area A make the same angle iJjq with respect to a± (see 
Fig.3(a)); on the other hand, when s = — 1, the angle 
increases proportionally to (see Fig.4(a)). 

In the calculation, we consider disclinations with two 
charges, s = d=l. The disclination area A and the mem- 
brane thickness h are set to be A = 4007T jam 2 and 
h = 4 /im by referring to the material constants of actual 
liquid crystal membranes [26]. The values of the elastic 
constants K\ and K3 are those of 4-methoxybenzylidene- 
4 / -butylaniline (MBBA), p-azoxyanisole (PAA), and 4- 
pentyl-4 / -cyanobiphenyl (5CB) presented in Table I [25]. 
The constant ipo in Eq.(3) is fixed to be zero without loss 
of generality, bearing in mind the fact that other choices 
of i/jq lead to no essential change in the conclusion of this 
work. 



III. RESULT 

Figure 2 shows plots of the disclination energy of 
MBBA as a function of the core position r c . The vertical 
axis shows the difference in Fd between the bump and 
the flat plane, i.e., Fd(a ^ 0) — Fd(a = 0). The discli- 
nation charge s is s = +1 in (a) and —1 in (b), and the 
parameter a ranges from a = to 1.2 in both (a) and 
(b). For s = +1, Fd takes the minimum at r c = ro for all 
values of a, indicating that a disclination is stable at the 
inflection point of the bump r = ro independently of the 
bump height. On the other hand, the plot for s — — 1 
exhibits a transition from a single downward peak struc- 
ture (a < 0.7) into a double- well structure (a > 0.7) such 



TABLE I: Elastic constants of MBBA, PAA, and 5CB re- 
ported in Ref.[25] 



Elastic constants 


MBBA 


PAA 


5CB 


(PN) 


(25°C) 


(122°C) 


(26°C) 


K x 


6 


6.9 


6.2 


K 3 


7.5 


11.9 


8.2 
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FIG. 2: Disclination energy in deformed MBBA films as a function of the core position r c . The parameter a characterizing 
the bump height ranges from to 1.2 (as indicated), (a) Energy of a positively charged disclination (s = +1). The magnitude 
of a downward peak at r c = ro monotonically increases with a. (b) Energy of a negatively charged disclination (s = — 1). 
Transition from a single- well to a double- well structure is observed at a = 0.7. 



that ro becomes unstable. The contrasting behaviors of 
Fd for s = +1 and —1 are attributed to surface curvature 
effects on the director configuration close to the core, as 
will be shown in Sec. IV. For PAA and 5CB, we have ob- 
tained similar profiles of Fd exhibiting a downward peak 
at r c = ro for s = +1 and a double- well structure for 
s = —1, whereas the peak heights differ slightly from 
those shown in Fig. 2 because of quantitative differences 
in K\ and K$. We have also confirmed that the peak 
positions (both downward and upward) are independent 
of our choice of r^o defined by Eq.(3). 

We emphasize that the energy scales of the potential 
wells shown in Figs. 3 (a) and 3(b) exceed the thermal en- 
ergy scale at ambient temperature. For instance, the 
depth of the single well for s = +1 and a = 1.0 given in 
Fig. 3(a) is approximately 5.7 x 10 _18 J; this is of the or- 
der of 10 3 k B T at T = 25°C. This fact suggests that the 
potential well at r = ro overwhelms the thermal excita- 
tion that causes a random movement of disclinations over 
the surface. As a result, an ensemble of trapped disclina- 
tions forms a halo with radius r c around the top of the 
Gaussian bump, provided that the disclination density is 
sufficiently low to ignore their repulsive interaction. The 
same scenario applies to disclinations with a charge of 
— 1, implying the formation of doubly consentric halos 
observed for a > 1.0. 



Figures 3(a) and 3(b) illustrate how the configuration 
varies with a shift in the core position r c from r c = 
(a) to r c = ro (b). In Fig. 3(a), n aligns radially around 
the top of the bump, exhibiting rotational symmetry such 
that n || a\ and n a,2 = at every point. This symmetry 
breaks down for r c > 0. By increasing r c , the perimeter 
of the disclination distorts slightly and becomes elliptical, 
and then, the relationship n • a2 = no longer holds (see 
Fig. 3(b)), except along the four specific line segments 
(6 = 0, 77, 7r, ■ 3 ^ L ). The most important observation in 
Fig. 3(b) is a suppression of the splay deformation in the 
shadowed region. The area of the corresponding region 
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IV. DISCUSSIONS 

The physical origin of a downward peak in the plot of 
Fd for s = +1 (see Fig. 2 (a)) is a curvature-induced shift 
in the director configuration around the disclination core. 



FIG. 3: Configurations of the director n around the core of 
a disclination with a charge of +1. With a shift in the core 
position from r c = (a) to r c = ro (b), the radial symmetry 
of the director configuration breaks down such that the splay 
deformation in the shadowed region in (b) is suppressed. 
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FIG. 4: Configurations of n in a disclination with a charge of —1. The core positions r c are r c = in (a), < r c < ro in (b), 
and r c = ro in (c). In the light- and dark-shadowed region, the splay deformation is suppressed and enhanced, respectively. 



is maximized at r c = ro , as a result of which a downward 
peak Fd arises at r c = ro at the expense of a slight splay 
enhancement within a limited region around 9 = 

Parallel arguments to the preceding paragraph account 
for the double- well structure of Fd for s = — 1 that is 
shown in Fig.2(b). Figures 4(a)-(c) show the variation 
of the n configuration under s = — 1 and a > 0.7, each 
figure corresponding to the core position: r c = (a), 
< r c < ro (b), and r c = ro (c). The director con- 
figuration in (a) exhibits a four-fold rotational symmetry 
that breaks down as r c increases. Symmetry breaking in- 
duces a suppression of the splay deformation in the light- 
shadowed regions shown in Fig. 4(b) and 4(c) in a manner 
similar to the case of s = +1. In contrast, the splay de- 
formation in the dark-shadowed regions is enhanced with 
an increase in r c . These two competing effects result in 
the double- well structure of Fd for s = — 1 cases. That 
is, Fd decreases with r c when the suppression effect ex- 



ceeds the enhancement, and it increases in the opposite 
situation. The degree of splay enhancement in the dark 
regions is maximized at r c = ro; this explains why ro is 
the most unstable point for s = — 1 disclinations. 

We have numerically confirmed that the splay term 
K\{D • n) 2 in Eq.(l) plays a dominant role in describing 
the r c - (and a-) dependences of Fd. Figures 5(a)-(b) 
show contour plots of the splay term K\(D • n) 2 within 
the disclination area of s = — 1. We observe that with an 
increase in r c , the local splay energy decreases at points 
near = tt/2 or 37r/2 but increases near = or 7r; this 
is in agreement with our previous discussions. A more 
detailed version of the study that covers the variations 
of the disclination size A and the initial angle ipo will be 
presented elsewhere. 
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FIG. 5: Contour plots of the splay term K\(D • n) 2 of discli- 
nations with s = — 1 placed at r c = (a) and at r c = ro 
(b). Light (dark) gray indicates a small (large) value of the 
term at the position. With an increase in r c , the splay term 
decreases around the line segments of = tt/2 and 3tt/2 but 
increases around those of = and ir. 



V. SUMMARY 



In conclusion, we investigated the elastic energy of 
disclinations that arise in a curved nematic membrane 
with a Gaussian bump geometry. When a < 0.7, the 
disclination energy Fd exhibits a downward peak at r = 
ro independently of the sign of the disclination charge, 
suggesting a curvature-driven halo formation of disclina- 
tions around the top of the bump. Furthermore, when 
a > 0.7, Fd for s = —1 exhibits a double- well structure 
in which the point r = ro becomes unstable in contrast 
to the case of s = +1; this implies doubly concentric ha- 
los peculiar to negatively charged disclinations. In both 
cases of s = ±1, the well depth overcomes the thermal 
energy scale at room temperature, suggesting the exper- 
imental feasibility of our theoretical predictions. 
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Appendix A: Metric tensor g^ v — ■ a v 

In this Appendix, we derive the metric tensor g^ v = 
a^-a v in terms of the covariant basis vectors a\ = dR/dl 
and a2 = dR/dO introduced in Sec. II. We start with the 
basis vectors e\ = dR(r,(j))/ dr and €2 = dR(r,(j)) / d(j) 
associated with the polar coordinate system (r, <p) on the 
flat plane. In terms of e\ and e2, the metric tensor g^ v = 
e M • e v becomes 



9nv 



f(r) 
r 2 



a 2 r 2 



> /(r) = l + -5-exp - 3 .(Al) 



We use a coordinate transformation from (r, <p) to (1,0) 
by shifting the coordinate origin from r = to r = r c . 
This is accomplished by multiplying both sides of g^ v by 
the Jacobian f3 defined by 

B _f dr/dl d<f>/dl\ ( 
£ ~ \ dr/dO d<f>/d6 J ' 

whose elements can be co mputed by using t he follow- 
ing relationships: r = yVj? + 2r c l cos 9 + I 2 and 4> = 
tan -1 {7 sin 9/(r c + lcos9)}. Finally we obtain 



9nv = a lx -a v = §g^(§) 1 

_ 1 ( f(r)p 2 + q 2 lpq(l " /W) 
" r 2 I lpq(l-f(r)) l\f(r)q 2 +p 2 ) ' 



(A3) 



where p = r c cos# + I and q = r c sinO. In particular, 
when a = (i.e., a flat plane), f(r) = 1, and thus, we 



1 

I 2 



(A4) 



as expected. 



Appendix B: Contravariant components n l 

We present the derivation of the contravariant compo- 
nent n % of the vector n = n % CLi in terms of a^. First, the 
relative angle between n and a\ is denoted as ip. Sim- 
ilarly, the angle between a\ and a 2 is denoted as It 
follows from Eq.(A3) that 



cos = 
sin^ 



a i • a 2 

|ai||a 2 | 
\ai x a 2 \ 



012 



/gH\/922 



\a>i\\a>2\ V^iiv^22 
where g = det^^]. We thus have 
1 



(Bl) 
(B2) 



cos (\£ — ijj) — 



V '011022 



(g 12 cos ip -\- *sj~g sin ip) , (B3) 



which leads us to the explicit forms of the covariant com- 
ponents rti = n • cti such as 



ri\ = |n||ai| costp = y/gu cos^, 

ri2 = \n\\a2\ cos (^ — ip) 
i 

1(012 cos ^ + ^/gs'mip). 



'9ii 



(B4) 
(B5) 



Finally, the formulas n J = n^* J with g^ ■ 
employed to obtain 

1 ( i 312 . ,x 
(cos w sinw , 

0ii Vg 
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